Long Time Algebraic Relaxation in Chaotic Billiards 
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The long time algebraic relaxation process in spatially periodic billiards with infinite horizon is 
shown to display a self-similar time asymptotic form. This form is identical for a class of such 
billiards, but can be different in an important special case. 
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Diffusion and transport of a population of particles 
whose motion occurs in an infinite periodic array of scat- 
ters is a basic problem of continuing interest (e.g., Refs. 
[jjj and 0]). Here we consider billiards, i.e., systems in 
which a point particle moves in straight line orbits with 
constant velocity, executing specular reflection from fixed 
boundaries (e.g, Refs. ||, ||, || ||). Several examples are 
shown in Fig. [I] and described in the figure caption. These 
examples all have channels in which a particle traveling 
in the proper direction (either exactly horizontally or ver- 
tically for the examples in Fig. never experiences re- 
flection; i.e., these are infinite horizon billiards. Particle 
motion without reflection occurs only for a zero mea- 
sure set of initial conditions, but it represents a possible 
source of deviation from classical diffusive behavior. De- 
viations from classical diffusive behavior (in particular, 
superdiffusive transport) are of great interest in a vari- 
ety of physical situations 0, [| [t|, ||, 0, 0, |lj|, |l3) and 
are often associated with the simultaneous presence of 
chaotic regions and invariant phase space surfaces [i.e., 
Kolmogorov-Arnold-Moser (KAM) surfaces]. In those 
cases the anomalous transport is associated with the 
"stickiness" jlj, [H| of the KAM surfaces; chaotic 
particles near these KAM surfaces tend to remain near 
them for long periods of time during which they experi- 
ence long flights Jjj . An analogous phenomenon is present 
in the infinite horizon billiards of Fig. [l] in that particles 
in the channels experience long flights if their direction 
of travel is nearly aligned with the channel. While KAM 
stickiness remains a rather difficult phenomenon to fully 
analyze, infinite horizon billiards, like those in Fig. |l|, are 
much simpler to study and offer possible insights into the 
general phenomenon of nonclassical diffusive transport. 

Our work on infinite horizon billiards will concentrate 
on the relaxation of particle distributions to their invari- 
ant long-time asymptotic form. In particular, we will 
study how a distribution, initially having no particles in 
phase space regions corresponding to long flights, repopu- 
lates these regions. We find, using analytical and numer- 
ical methods, that this repopulation process occurs in a 
self-similar manner That is, at long time the relaxing 
distribution function assumes a particular invariant form 
when expressed in terms of a properly scaled variable. 
Furthermore, this distribution function is the same for 
all the billiards in Figs. |l](a)-(d), but is different for the 
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FIG. 1: Infinite horizon billiards, (a) Circular scatterers 
of radius R whose centers are located on a square grid of 
periodicity length L > 2R. This billiard is also referred to 
as the Sinai billiard ||. (b) The circular scatterers of (a) are 
each given a random displacement A that is uniform in the 
disc |A| < Ao < (L — 2R)/2. (c) Randomly oriented squares 
of edge length d where the centers of the squares are located 
on a square grid of periodicity length L > y/2d. (d) A channel 
whose walls are circular arcs each subtending the same angle 
which is less than n radians, (e) Similar to (d) but the circular 
arcs subtend n radians (i.e., they are semicircles); orbits in 
this billiard can be mapped to orbits in the stadium billiard 
of Bunimovich || by reflection at the horizontal dashed lines 
in the figure. 



billiard in Fig. |(e). The billiard in Fig. ||(e) is of partic- 
ular interest because its infinite flight invariant set is, in 
an appropriate sense, more sticky than are the invariant 
sets in the other examples. 

It is important to note that the billiards in Fig. |l| 
all admit an invariant particle probability distribution 
function (PDF). In particular, considering a monocner- 
getic ensemble of particles, a particle distribution func- 
tion that is uniform in the accessible area of the billiard 
and isotropic in the angle 9 of the particle velocity vector 
v is time invariant. (Here 6 is defined by \v-y \ — \v\ cos 8, 
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< 6 < tt/2 and y is a unit vector in the vertical di- 
rection.) Our main concern in this paper is how this 
invariant distribution, isotropic in 9, is approached. 

Since the channels where infinite flights occur are of 
most interest, we are primarily concerned with the man- 
ner in which particles enter and leave the vicinity of the 
vertical infinite flight orbits. For the purpose of exposing 
the essential features of this problem, we consider the re- 
laxation of an initial distribution, P(x,y,9,0), that has 
no particles in the vicinity of the infinite flight orbits. 
In particular, we take P{x, y, 9, 0) to be zero outside one 
of the size L cells indicated in Fig. E], while within the 
accessible part of such a cell 
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where if is a normalization constant, 9 max <C 1. 
Our main result is that for large time P(9, t) = 
J P(x,y,9,t)dxdy, where the integral is over (x,y) re- 
stricted to the vertical channels, approaches a scaling 
form depending on the single variable </> = v9t/W. That 
is, for large t and </> bounded, 



(2) 



Furthermore, the form of the scaling function on the right 
hand side of (||) is the same for the billiards in Figs. |l|(a)- 
(d). We call this function So(4>). For the billiard of 
Fig. |l|(e), however, the scaling function, denoted Si (</>), 
takes another form. In the definition <fi = v9t/W, W is 
defined as the region of free vertical flights [for billiard 
(a) W = L- 2R, for billiard (b) W = L- 2{R+ A ), for 
billiard (c) W is the distance between the dashed circles 
(W = L — V^d), and for billiards (d) and (e) W is as 
shown in Fig. [|.] 

Figure ||(a) shows histogram approximations to P(9 7 i) 
plotted versus <fi = v9t/W for the billiard of Fig. [l](a) 
with L — 3R at four times (+, x, *, and □ correspond 
to t = 9W/v, 27W/v, 8lW/v, and 243W/v respectively). 
These data are obtained by computing the orbits of a 
large number of initial conditions chosen randomly from 
the distribution function (Q). Results similar to those in 
Fig. ||(a) are also obtained for the billiards of Figs. |l|(b, c, 
and d). In all these cases, at long time, the distribution 
increases linearly with <fi from </> = 0, until some critical 
value (j)Q past which the distribution is constant in </>, 



P(M) = S o (0) 
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In contrast to the results [e.g., Fig. ||(a)] for the billiards 
of Figs. |(a-d), the billiard of Fig. |(e) gives a very differ- 
ent time asymptotic distribution. This is shown by the 
histogram approximations (Fig. ||(b)) to P(9,t) obtained 
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FIG. 2: (a) The time evolution of an initially uniform angu- 
lar distribution of particles in the 9 interval [.3,7r/2] for the 
billiard of Fig. [I](a) with L/R = 3. The distribution is shown 
at four times; + at t = 9W/v, x at 27W/v, * at 81W/v, and 
□ at 243W/v. (b) Distribution function for the billiard of 
Fig. |l](e) with W/ R = 1 at the same four times as in (a) and 
with the same initial 9 distribution as in (a). The solid line 
in (b) is obtained from the theory, Eqs. (^) and (|^). 



for the billiard of Fig. |l|(e) with W/R = 1 and the same 
four times as in Fig. ||(a) . Note that in this case the long 
time distribution function is identically zero for (j) < 1/2. 

We now show how the time asymptotic distributions 
in Figs. § arise. We start with the distribution Eq. (||) 
[Fig. ||(a)] applicable to the billiards in Figs. |l|(a-d). As 
we will subsequently see, for all the billiards in Figs. |l|(a- 
d), the scattering of a long flight upon reflection from a 
channel wall leads to a much more drastic change in the 
angle of a particle's velocity vector than in the case of 
the billiard in Fig. 0(e). 

For example, in the case of the billiard in Fig. |l|(a), the 
angular deflection for small 9 is typically of order 9 1 ^ 2 , 
which is much larger than 9. In the case of the billiards in 
Fig. [l](b-d), a particle moving nearly parallel to a channel 
axis is scattered by an angle of order one. Furthermore, 
after a large deflection, the orientation of the particle's 
velocity vector is rapidly randomized by a succession of 
many reflections which, since the particle is no longer in 
a long flight, occur in a relatively short time. These con- 
siderations lead us to a model for the cases in Figs. |l|(a)- 
|l|(d) in which we adopt the model hypothesis that, when 
a particle in a long flight suffers a collision with a billiard 
wall, the orientation of its velocity vector is randomly 
scattered with uniform probability density in [0, 27r]. We 
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FIG. 4: (a) A one-bounce encounter with the wall resulting 
in the reversal of the particle's vertical velocity component, 
(b) A two-bounce encounter with the wall resulting in an un- 
changed direction of the particle's vertical velocity compo- 
nent, (c) M(e) verses e for e ma x = 2i?tan#' <C 1. 



FIG. 3: The evolution of P c (0,t) and P(0,t). 



wish to determine the evolution from the initial condi- 
tion (0) in the case (9 max *C 1. This initial distribution 
is equal to the initial distribution Pb(x,y,9,0) = K for 
all 9 minus the initial condition P c (x,y, 9, 0) = K for 
< 9'max and otherwise. The distribution p, remains 
unchanged when it is evolved forward in time (it is an 
invariant distribution). Thus, to find the evolution from 
initial condition (0), we can determine the evolution from 
P c , and then subtract it from K. The long time evolu- 
tion from P c can be found by considering the time at 
which particles are scattered. Consider, for example, the 
billiard of Fig. 0(d), and a particle with a small initial 
o . Suppose the particle is located in the channel at 
a distance Ax from the boundary of the channel with 
which it will collide [left or right vertical dashed line in 
Fig. 0(d)]. If Ax < vtsm9 Q = vt9 D , the particle scat- 
ters before time t; if Ax > vt9 a , it does not scatter. For 
particles in the channel, Ax < W, so every particle with 
8 > W/(vt) must have scattered at least once. We as- 
sume that t > W/(v9 max ) = t a . Since (9 max is small, 
the scattered particles contribute a small positive value 
of order # max to P c (9,t) in 9. Thus P c (9,t) is small (i.e., 
of order max ) for 9 > W/vt. For 9 a < Ax/(vt), t > t D , 
the particle has not yet scattered. Assuming that the 
initial spatial distribution of particles in the channel is 
uniform, the fraction of particles with initial angle 9 Q 
that have scattered is 9 vt/W. Thus 



flights, Fig. 0(b). When an orbit in a long flight en- 
counters the channel wall it may experience either one 
reflection (Fig. 0(a)) or two reflections (Fig. 0(b)). A ge- 
ometrical analysis || shows that for 9' < 1, the relation 
between the angle 9' before the encounter and the angle 9 
after the encounter is 9 = M(e)9' , where e is the distance 
illustrated in Figs. 0(a,b) and M(e) is the piecewise lin- 
ear function shown in Fig. 0(c). Note from Fig. 0(c) that 
the trajectory angle can increase or decrease at most by 
a factor of 3 upon an encounter with the wall. Thus a 
long flight is again followed by a long flight, in contrast 
to the billiards of Figs. 0(a-d). That is, the infinite flight 
invariant set for the billiard of Fig. 0(e) is more sticky 
than are those of the other billiards, Figs. 0(a-d). Note 
also that, for 9 <C 1, e at the next encounter with the 
end of the channel is very sensitive to a small percent- 
age change in 9. This fact, together with the form of 
M(e) shown in Fig. 0(c), leads us to adopt a model for 
our subsequent analysis whereby upon an encounter with 
the wall 8' — ► 9 = M9' where M is chosen randomly with 
uniform probability density in the interval (1/3,3). As 
previously noted S\{<j>) = for <p < 1/2. This can be 
understood as follows. P(9,0) = for 9 < 9o, and on 
subsequent wall encounters 9 can decrease at most by a 
factor of 3. Therefore, after n encounters the minimum 
angle is 9 n = 9$/3 n and the minimum time for these n 
encounters is t n = YZ~=o W K v6 n<) ~ (W/v)3 n /(29 a ) for 
n 1, yielding a minimum value of 1/2 for <f>. 



P c (0,i) = 



K(l - 9vt/W) for 9 < W/(vt), 
for 9 > W/(vt), 



(4) 



where we have neglected the small, order 9 ma x, contri- 
bution to P c (8,t) from scattered particles. Subtracting 
P c from Pf, as illustrated in Fig. ||, we obtain the time 
asymptotic form in Eq. (0). 

We now consider the evolution of initial condition (0) 
in the billiard of Fig. 0(e) for orbits experiencing long 



Let Q(9,t)d9 be the rate at which particles are scat- 
tered into the channel with angle between 9 and 9 + d6. 
Particles that enter the channel with small angle 9 re- 
main in the channel for the time W/(vsm0) = W/{v9). 
Then for large t and small 9, since the number of parti- 
cles in the interval 9 to 9 + d# at time t is approximately 
St(4>)d9, we have S^) = J*_ w/{v9) Q(9,t')dt' . Since 
the left hand side is only a function of <fi, Q(9,t) must 
have the asymptotic scaling form Q(9, t) — (v6/W)C((f>). 
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Thus 



(5) 



To obtain an equation for C{4>) consider a particle in a 
long flight that, at time t' , has just suffered an encounter 
with a channel wall. Denote its angle as a result of this 
encounter by 9' <§C 1 (see Figs. |](a,b)). The particle then 
experiences a long flight of duration W/{v9') before it 
next encounters a channel wall. According to our model, 
the velocity angle 9 just after this encounter at time t is 
9 = MO' where M is random with uniform probability 
density in the interval (1/3,3). During a time interval t 
to t+dt the number of particles scattered into the interval 
9 to 9 + d9 is Q{9, t)A9dt = 3/8 f* /3 [Q(9', t')66'At'] AM, 

where 6' = 6/M, t' =t - $r, and 3/8 = 1/(3 - 1/3) is 
the probability density of the random variable M. Now 
substituting Q(6, t) — (v6/W)C(<f>) and making a change 
of the integration variable m = 1/M, we arrive at the 
following integral equation for C(<fi), 



C(4) = 3/8 C{m<t>-l)Am = L[C{(t>)] (6) 

A/3 

whose solution along with (||) determines Si. We nu- 
merically solve (||) by iteration of d+i(tp) = C[d((j))} 
starting with C (<f>) — if < 1 and C (4>) = 1 if <f> > 1. 
We use 10 4 grid points uniformly spaced in the 4> interval 
[0, 50] with Ci(4>) set to one for <p > 50 (this corresponds 
to the boundary condition C(cj>) — > const, for <f> — * oo). 
The results are virtually unchanged if the <fi interval is 
increased to [0, 100] or if the number of grid points is in- 
creased, and virtually identical results (i.e., convergence) 
were found for several i > 12. A result obtained from 
such a calculation is shown as the solid curve in Fig. || 
(i = 15) and is in good agreement with our numerical 
histograms computed from many particle orbits fl7|] . 

In conclusion, we have shown that the long-time repop- 
ulation of long flights in infinite horizon billiards proceeds 
in a self-similar manner [Eq. (0)]. The corresponding re- 
laxation of the distribution [which follows from (g)] is 
algebraic rather than exponential; e.g., the 9 integral of 
\P(6, t) — P(9, oo)| decays like l/t. Furthermore, the self- 
similar forms found for a class of billiards that include 
those of Figs. |l|(a-d) are the same, Eq. (||). On the other 



hand, the special case of the billiard of Fig. |T](e) is shown 
to yield a different time-asymptotic self-similar form, and 
this may be ascribed to the greater stickiness of the in- 
variant set of this billiard. 
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